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This is a colloquium-style introduction to two electronic processes in a carbon monolayer 
(graphene), each having an analogue in relativistic quantum mechanics. Both processes cou- 
ple electron-like and hole-like states, through the action of either a superconducting pair potential 
or an electrostatic potential. The first process, Andreev reflection, is the electron-to-hole conver- 
sion at the interface with a superconductor. The second process, Klein tunneling, is the tunneling 
through a p-n junction. The absence of backscattering, characteristic of massless Dirac fermions, 
implies that both processes happen with unit efficiency at normal incidence. Away from normal 
incidence, retro-refiection in the first process corresponds to negative refraction in the second pro- 
cess. In the quantum Hall effect, both Andreev reflection and Klein tunneling induce the same 
dependence of the two-terminal conductance plateau on the valley isospin of the carriers. Existing 
and proposed experiments on Josephson junctions and bipolar junctions in graphene are discussed 
from a unifled perspective. 
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I. INTRODUCTION 



In October 1920 the inventor of special relativity, Al- 
bert Einstein, traveled to Leiden to meet with the discov- 
erer of superconductivity, Heike Kamerlingh Onnes. A 
photograph of a blackboard (Fig. [T]) records one of their 
discussions. The two physicists had much to discuss, but 




FIG. 1 Albert Einstein, Paul Ehrenfest, Paul Langevin, Heike 
Kamerlingh Onnes, and Pierre Weiss at a workshop in Leiden 
(October, 1920). The blackboard discussion, on the Ha l l effec t 
in super conductors , has b een reconstructed by ISauerl l|2007f ). 
See also IVan DelftI i 20061 ) for the historical context of this 
meeting. 



they would have found little common ground in the two 
topics closest to their hearts, since superconductivity is 
essentially a nonrelativistic phenomenon. 

Relativistic contri butions to the superconducting pair 
potential, studied bv lCapelle and GrossI (|l995l[l999l ). are 
a small correction of order (vp/c)'^ [Fermi velocity over 
speed of light squared] . Fully relativistic phenomena such 
as particle-to-antiparticle conversion by a superconduc- 
tor have remained pure fiction. Some of this fiction is 
now becoming science in a material first isolated a few 
years ago b y Andre Geim and his group at Manchester 
University i Novoselov et al.l . l2004) . 




The material, called graphene, is a mono-atomic layer 
of carbon atoms arranged on a honeycomb lattice. Upon 
doping, electrons and holes move through the layer with 
a velocity v = lO^m/s which is only a small frac- 
tion of the speed of light. And yet, this velocity is 
energy independent — as if the electrons and holes 
were massless particles and antiparticles moving at the 
speed of light. As demonstrate d in transport rneasure - 
ments bv lNovoselov et aL ] (l2005^ and lZhang et al. |(|200i, 
and in spectroscopic measurements bvl Zhou et al.l (I2OO6D 
and Bostwick ct al. ( 2007), the electronic properties of 
graphene are described by an equation (the Dirac equa- 
tion) of relativistic quantum mechanics, even though the 
microscopic Hamiltonian of the carbon atoms is nonrel- 
ativistic. While graphene itself is not superconducting, 
it acquires superconducting properties by proximity to a 
superconductor. We therefore have the unique possibility 
to bridge the gap between relativity and superconductiv- 
ity in a real material. 

For example. Fig. [5] shows two superconducting elec- 
trodes on top of a carbon monolayer. T he super c urrent 
measured through this device bv .Heersche et al.l ( 20071 ) 
is carried by massless electrons and holes, converted 
into each other by the superconducting pair potential. 
This conversion process, known as Andreev reflection 
( Andreevl . ll964 ). is describ ed by a supercoii ducting vari- 
ant of the Dirac equation ( Beenakkerl l200 d) . 

In this Colloquium we review the unusual physics 
of Andreev reflection in graphene. For a broader 
perspective, we compare and contrast this coupling 
of electrons and holes by a superconducting pair po- 
tential with the coupling of electrons and holes by 
an electrostatic potential. The latter phenomenon 
is called Klein tunneling (ICheianov and Fal I2OO6I: 
iKatsnelson. Novoselov. and Geiml . 20061) with reference 
to relativistic quantum mechanics, where it represents 




FIG. 3 Band structure E{kx, ky) of a carbon monolayer. The 
hexagonal first Brillouin zone is indicated. The conduction 
band [E > 0) and the valence band {E < 0) form conically 
shaped valleys that touch at the six corners of the Brillouin 
zone (called conical points, or Dirac points, or if-points). The 
three corners marked by a white dot are connected by recipro- 
cal lattice vectors, so they are equivalent. Likewise, the three 
corners marked by a black dot are equivalent. In undoped 
graphene the Fermi level passes through the Dirac points. 
(Illustration by C. Jozsa and B. J. van Wees.) 



the tu nneling of a particle into the Dirac sea of antipar- 
ticles (lKleiniril929D . Klein tunneling in graphene is the 



tunneling of an electron from the conduction band into 
hole states from the valence band (which plays the role 
of the Dirac sea). 

The two phenomena, Andreev reflection and Klein 
tunneling, are introduced in Sees. IIIII and llVi re- 
spectively, and then compared in Sec. |Vl But 
first we summarize, in Sec. |lTl the special prop- 
erties of graphene that govern these two phenom- 
ena. More comprehensive reviews of graphene have 



been written by Castro Neto. Guinea, and Pered ( 


2006^; 
2007i): 


Castro Neto et al. (l2007l): iGeim and Novoselovl ( 


Gusvnin. SharaDOv. and Carbotte 


(I2OO7I); KatsnelsonI 


(2007): Katsnelson and Novoselov 


(120071). 



II. BASIC PHYSICS OF GRAPHENE 

A. Dirac equation 

The unusual band structure of a single layer of 
graphite, shown in Fig. [31 has been known for 60 years 
(| Wallace! . fl947l) . Near each corner of the hexagonal first 
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FIG. 4 Honeycomb lattice of a carbon monolayer. The unit 
cell contains two atoms, labeled A and B, each of which 
generates a triangular sublattice (open and closed circles). 
The lattice constant a is \/3 times larger than the carbon- 
carbon separation of 0.142 nm. The reciprocal lattice vector 
K has length 4n/3a. The edge of the lattice may have the 
armchair configuration (containing an equal number of atoms 
from each sublattice) , or the zigzag configuration (containing 
atoms from one sublattice only). Dashed circles and bonds 
indicate missing atoms and dangling bonds, respectively. The 
separation W of opposite edges is measured from one row of 
missing atoms to the opposite row, as indicated in the figure. 



Brillouin zone the energy E has a conical dependence on 
the two-dimensional wave vector k = (k^, ky). Denoting 
hy 6k — k — K the displacement from the corner at wave 
vector K, one has for \Sk\a <C 1 the dispersion relation 



\E\ = hv\6k\. 



(1) 



The velocity v = ^^/Sra/U « lO^m/s is proportional 
to the lattice constant a — 0.246 nm and to the nearest- 
neighbor hopping energy r w 3 eV on the honeycomb 
lattice of carbon atoms (shown in Fig.|3|). 

The linear dispersion relation ^ implies an energy- 
independent group velocity Wgioup 

EE \dE/ndk\ ^ V oi 
low-energy excitations {\E\ ^ r). These electron ex- 
citations (filled states in the conduction band) or hole 
excitations (empty states in th e valence b and) there- 
fore have zer o effect i ve m ass. ISemeno^ (|l984l ). and 
lOiVincenzo and Meld (|1984| ) noticed that — even though 
?; <C c — such massless excitations are governed by a 
wave equation, the Dirac equation, of relativistic quan- 
tum mechanics: 



ihv ( ^ " 

dx + idy 



= E 



(2) 



[The derivat ion of this equa tion for a carbon monolayer 
goes back to lMcClurd (|1956[) .] 

The two components "^^a and ^'b give the amplitude 
^^(r)e'''^ '' and ^B(r)e'''^ '' of the wave function on the 
A and B sublattices of the honeycomb lattice (see Fig. [4]). 
The differential operator couples ^'^ to but not to 
itself, in view of the fact that nearest-neighbor hopping 
on the honeycomb lattice couples only ^-sites with B- 



^ tpB, Ei)B=T ^ ■0/ 



(3) 



ehbo 



E^hbo 



In a more concise notation, Eq. ^ may be written as 



Vp • (Tip — Eip^ 



(4) 



with p = —ih{d/dx, d/dy) the momentum operator in 
the x-y plane and cr = {(Tx, (Jy, az) the vector of Pauli 
matrices acting on the spinor -0 = {'^a, ^b)- (For later 
use we define ctq as the 2x2 unit matrix.) The spin 
degree of freedom described by the Pauli matrices at is 
called the "pseudospin" , to distinguish it from the real 
electron spin. 

This two-dimensional Dirac equation describes states 
with wave vector k in the valley centered at the cor- 
ner of the Brillouin zone with wave vector K = 
{4TT/3a)x. The valley at the opposite corner at —K 
produces an independent set of states with amplitudes 
^'ji{r)e-'^ '^ and *^(r)e-*^ '' on the A and B sub- 
lattices. The two components ^'j^ and '^'g satisfy the 
same Dirac equation (g]) with px ~Px- The spinor 
^ = (VI^A, -^'s, containing both valleys there- 
fore satisfies the four-dimensional Dirac equation^ 



vp ■ cr 







vp ■ cr 



* = E^. 



(5) 



This differential equation represents the low-energy and 
long-wave length limit of the difference equation ([3|) in 
the tight-binding model of graphene. 

For a compact notation, we make use of a second set 
of Pauli matrices t = {tx, Ty,Tz), with tq the 2x2 unit 
matrix, acting on the valley degree of freedom (while cr 
and (Jo act on the sublattice degree of freedom). Eq. ([5]) 
may then be written as 



H{A)^ = E^, 

H{A) = v[{p + eA) • cr] (g) To 



Uao (g) To, 



(6a) 
(6b) 



where for generality we have also included external elec- 
tromagnetic fields (with scalar potential U and vector 
potential A). Electromagnetic fields do not couple the 
two valleys, provided that the fields vary smoothly on 
the scale of the lattice constant. 

To conclude this subsection we briefly comment on 
the quantum-relativistic analogue of Eq. ([5]), referring 



Next-nearest-neighbor hopping contributes second-order spatial 
derivatives, which arc of higher order in a\Sk\ and may therefore 
be neglected in first approximation. 

The valley-isotropic representation ((5J of the four-dimensional 
Dirac equation (with two identical 2x2 subblocks) is used to 
write boundary conditions in a compact form (see Sec. III.C^ . 
Other representations (with two unequal subblocks) are common 
in the literature as well, and one should be aware of this when 
comparing formulas from different papers. 
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to iGusvnin. Sharapov. and Carbottd ( 20071) for a more 
extensive discussion. In three dimensions, and with a 
change of sign of one of the two subblocks vp - (t, Eq. ([5|) 
represents the Dirac (or Dirac-Weyl) equation of massless 
neutrinos, with v the speed of Ught. The vaUey degree 
of freedom corresponds to the chirahty of the neutrinos, 
which have left-handed or right-handed circular polariza- 
tion (corresponding to the opposite sign of the two sub- 
blocks). In two dimensions the relative sign of the two 
subblocks can be changed by a unitary transformation, 
so the distinction between left- or right-handedness can- 
not be made. Electrons in graphene are called "chiral" 
because their direction of motion is tied to the direction 
of the pseudospin. Indeed, the current operator 



vcr §5 To 



(7) 



is proportional to the pseudospin operator cr, so that an 
electron moving in the x or y-direction has a pseudospin 
pointing in the x or j/-direction. But because the pseu- 
dospin is two-dimensional, there is no analogue of circu- 
lar polarization and therefore there is no left- or right- 
handedness in graphene. 



B. Time reversal symmetry 

The time reverse of the state ^xe^^"^ + 
on the X = A,B sublattice is the complex conjugate 



-iK-r 



verse of the spinor ^ 
therefore has the form 



This implies that the time re- 
The time reversal operator T 



r 



/GO IN 
0-10 
0-100 

\1 Oy 



-{Ty (g) (Ty)C, 



(8) 



with C the operator of complex conjugation. Notice 
that the time reversal ope ration interchanges the valleys 
(|Suzuura and Andol . l200^ . 

The time reverse of the Dirac Hamiltonian (161) is 



TH{A)r-^ = H{-A) 



(9) 



As it should be, time reversal symmetry is preserved in 
the absence of a magnetic field. 

The Dirac Hamiltonian satisfies another anti-unitary 
symmetry, SH{A)S^^ = H{-A), with 



is called a symplectic symmetry, while the T-symmetry 
is called an orthogonal symmetry.'^ 

Because it is not the true time reversal symmetry oper- 
ator on the honeycomb lattice, the symplectic symmetry 
can be broken even in the ab sence of a magnetic field 
( Berrv and Mondrago Consider the following 

two perturbations SH of the Dirac Hamiltonian: 

• A mass term 6H — fi{r)az (X" t^, generated for ex- 
ample by a sublattice depe ndent potential in the 
substrate (|Zhou et al.l . l2007t ). 



• A valley-dependent vector potential, SH = 
ev\A( r) ■ tr] (g) T^, produced by straining the mono- 
layer (jMorozov et al.. . ,200&,Morpurgo and Guineal . 
l200fih .4 

In both cases, T6HT^^ — 6H, so time reversal symmetry 
is preserved, while S6HS^^ = —SH, so the symplectic 
symmetry is broken. 

Whether it is the T-symmetry or the 5-symmetry that 
governs a transport property depends on whether the 
scattering processes couple valleys or not. A smoothly 
varying electrostatic potential does not cause interval- 
ley scattering, so it is the presence or absence of the 
symplectic symmetry S that matters in this case. [For 
example, breaking of S destroys th e weak antilocaliza- 
tion e ffect, even if 7 " is pr e served (lAleiner and Efetovl . 
120061: iMcCann et all . 120061 : ISuzuura and Andd . 12002^ ] 
Andree v reflection at a su perconductor does couple the 
valleys ( Beenakkeil . 120061 ). so there it is the true time 
reversal symmetry T that matters. [For example, break- 
ing of T s uppresses the supercu rrent while breaking of S 
does not (jHeersche et all . l2007t) .] 



C. Boundary conditions 

The Dirac equation needs to be supplemented by a 
boundary condition o f the form ^ = My at th e edge 
of the graphene sheet ( McCann and Far"koll2004l ). Since 
edges are typically abrupt on the atomic scale, the bound- 
ary condition couples the valleys. Ignoring a possible 
local magnetization, we may assume that M commutes 
with T — meaning that the boundary condition itself 
does not break time reversal symmetry. The boundary 
condit ion then has the form ( Akhmerov and Beenakken . 
l2007bf ) 



^' = Af*, M ^ {u ■ t) (g) {n ■ cr). 



(11) 



1 ON 

-10 

1 

-1 Oy 



C = i{TQ g) (Jy)C. 



(10) 



This operator S does not interchange the valleys, unlike 
T, but like T it does invert the sign of p and a. The op- 
erator S, therefore, acts like a time reversal operator in a 
single valley. The 5-symmetry of the Dirac Hamiltonian 



A symplectic symmetry operator is an anti-unitary operator 
which squares to —1, while an orthogonal symmetry operator 
is an anti-unitary operator which squares to +1. Both T and 
S are anti-unitary (product of a unitary operator and complex 
conjugation), but = 1 while = —1. 

A ripple of diameter R and height H corresponds to a fictitious 
magnetic field of order \B\ ~ {Ti/ ea)H^ / B? , of opposite sign in 
the two valleys. 



FIG. 5 Location of the valley isospin v on the Bloch sphere 
for a zigzag edge (blue arrows) and for an armchair edge (red 
arrows). The solid and dashed arrows correspond to opposite 
edges. 



parameterized by a pair of three-dimensional unit vectors 
V and n. The vector n is constrained by n • ub = 0, to 
ensure that no current leaks out through the boundary 
(with normal n^, pointing outward). 

We give three examples of boundar ies 
(jBerrv and Mondragonl llQSTt iBrev and Fert i3. l2006aD : 

• A zigzag edge has either a = ^'a = or = 
^''g — 0, depending on whether the row of miss- 
ing atoms at the edge is on the A or B sublatticc 
(see Fig. U). The corresponding boundary condi- 
tion matrix M has v = zLz, n = z. Because op- 
posite zigzag edges lie on different sublattices, the 
angle $ between the vectors v on opposite edges 
equals tt — irrespective of the separation of the 
edges. 



^'^e-'^ '' = for 



• An armchair edge has '^xe^'^ 
X = A,B — so that the wave function vanishes on 
both sublattices. This requires i'-z = 0, n = zx 
ub- The angle $ = |iir|M^ -I- tt now depends on the 
separation W (as defined in Fig. [4]): $ = tt if 2W/a 
is a multiple of three and <i> ~ ±7r/3 otherwise. 

• Confinement by an infinite mass has u — z, n = 

z X Ub- 

The two eigenstatcs | + v) and | — i'') of i^* • t (defined by 
u-tIzLv) = ±|±i/)) are states of definite valley polariza- 
tion — parallel or antiparallel to the unit vector u. This 
vector is called the valley isospin, because it transforms 
under rotations in the same way as the real electron spin. 
It can be represented by a point on the Bloch sphere, see 
Fig.O When v points in the z direction, the polarization 
is such that the eigenstate lies entirely within one single 
valley. This is the case for the zigzag edge or for the in- 
finite mass confinement. When u lies in the x-y plane, 
the eigenstate is a coherent equal-weight superposition of 
the two valleys. This is the case for the armchair edge. 

The direction of u in the boundary matrix M plays a 
key role in a strong magnetic field, by selecting the valley 



FIG. 6 These two graphene flakes both have the same zigzag 
boundary condition: = ±rz ® a^'i'- The sign switches be- 
tween -|- and — at the armchair orientation (when the tangent 
to the bounda ry has an angle with the y-ax\s wh ich is a mul- 
tiple of 60°). (jAkhmerov and Beenakkeil . l2007cD 



polarization of edge states ( Akhmerov and Beenakke^ 
[200 7b) . Edge stat es in the lowest Landau level 
are valley polarized ( Abanin. Lee, and Levitovl . l2006t 
iBrev and Ferti g|. l2006bl ). but the Hall conductance is in- 
sensitive to the direction v of the valley isospin.^ In Sec. 
IV. CI we will see how Andreev reflection and Klein tun- 
neling both provide a way to measure the valley isospin 
in the quantum Hall effect. 

We conclude this discussion of boundary conditions 
with the constraint imposed by electron-hole symmetry. 
In the absence of an electrostatic potential {U = 0), the 
Dirac Hamiltonian ^ anticommutes with (8) cr^ . In an 
unbounded system, this implies electron-hole symmetry 
of the spectrum. (If ^E" is an eigenstate with eigenvalue 
E, then Tz (8) (Jz^ is an eigenstate with eigenvalue —E.) 
The electron-hole symmetry exists already at the level of 
the tight-binding model ^ [E ^ -E if *b 
so it is preserved by any boundary that is simply a ter- 
mination of the lattice (zero edge potential).^ The re- 
quirement that the boundary matrix M in Eq. (jlip com- 
mutes with Tz ® <Jz (needed to preserve the electron- 
hole symmetry) restricts M to either the zigzag form or 
the armchair form. As illustrated in Fig. [51 the zigzag 
form is th e rule while the a r mchair form is the exception 
(lAkh merov and Beenak"km2007cl) . 



D. Pseudo-diffusive dynamics 

Electrical conduction through a graphene sheet has 
unusual features when the Fermi level is at the Dirac 



^ The Hall conductance Gjj = 9^ 1^ is determined by the de- 
generacy factor g of the edge states. The celebrated "half- 
int ege 



Hall conductance Gh 



{n-\- 1/2) X ie^ /h measured 
by iNovoselov et al.l |200i) and IZhang et all ||2005| ') tells us that 
the lowest (n = 0) Landau level has spin degeneracy but no val- 
ley degeneracy (g = 2 rather than g = A). The direction of the 
valley polarization does not enter in Gh- 
' One mechanism that may produce an edge potential at 
a zigzag boundary (antiferromagnetic spin ordering) has 
been discussed in connectio n with graphene nanoribbons by 
ISon. Cohen, and Louil ||2006|1 . 
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FIG. 7 Electrostatic potential profile (red solid line) produc- 
ing two heavily doped graphene regions at the left and right 
and a weakly doped region (length L) at the center. The cen- 
tral region is undoped when the Fermi energy (blue dashed 
line) coincides with the energy of the Dirac point. Electri- 
cal conduction then proceeds via evanescent (= exponentially 
decaying) modes. 




valence band 



FIG. 8 Electron and hole excitations in the conical band 
structure of graphene (filled and empty circles at energies 
Ef ±e), converted into each other by Andreev reflection at a 
superconductor. 



point. Because of the vanishing density of states the 
transmission through a strip of undoped graphene (width 
W, length L in the current direction) occurs entirely 
via evanescent {— exponentially decaying) modes. For 
a short and wide strip there is a large number W/ L 3> 1 
of evanescent modes with transmission probability of or- 
der unity. In a remarkable coincidence/ the transmis- 
sion probabilities of the evanescent modes are the same 
as those of diffusive modes i n a disordered piece of m etal 
with the same conductance ( Tworzvdlo et al.1 . l2006l ). We 
will return to this "pseudo-diffusive" dynamics in Sec. 
IIII.Dl when we describe how supercurrent flows through 
undoped ballistic graphene in the same way as it does 
through a disordered metal. 

In preparation of that discussion, we examine here 
in a bit more detail the transmi ssion of evanescent 
modes through undope d graphene ( Katsnels"onl l2006ht 
iTworzvdIo et al.l . I2006D . Because the wave length at 
the Dirac point is infinitely long, the detailed shape of 
the electrostatic potential profile at the interface be- 
tween the metal contacts and the graphene sheet is not 
very important. We model it by the rectangular poten- 
tial shown in Fig. [T] The contact area is modeled by 
he avily do ped graphene [for more microscopic models. 



see 



Blanter and Martin (20071): [Robinson and Schomerud 



(|2007l ): ISchomerusI (|2007i )]. The Fermi level in Fig. [7]hes 
in the conduction band in the contact areas at the left 
and right and in the valence band in the central region. 
Conduction in this situation occurs via interband (Klein) 
tunneling, from conduction band to valence band, and wc 
will have much more to say about that in Sec. IIVI 

The special situation we are interested in here is when 
the Fermi energy coincides with the energy of the Dirac 



We say "coincidence" because we have no intuitive explanation 
for this correspondence. 



point in the central region. At that energy interband tun- 
neling goes over into intraband tunneling. For W/L 3> 1 
we do not need to know the individual transmission prob- 
abilities of the evanescent modes (which will depend on 
the boundary condition at y = 0,11^), but it suffices to 
know how many modes p(T)dT (counting all degeneracy 
factors) there are with transmission probabilities in the 
interval (T, T + dT). The result is 



P{T) 



(12) 



with g the conductance in units of e^/h. We call 
the dynamics pseudo-diffusive because the distribution 
( happens to coincide with the known distribution 



DorokhovL 11984 for diffusion modes in a disordered 



metal having the same dimensionless conductance \. 

The value for q in Eg. (1121) has been confirmed exper- 
imentally bv lMiao"et al.. (20071) . To test for the bimodal 
shape of the distribution p{T) one would need to measure 
the shot noise at the Dirac point. The Fano factor (ra- 
tio of shot noise power and mean current) should equal 



(jTworzydlo et al 



power I 

Illooi 



F =1- 



jlT^p{T)dT ^1 
SlTp{T)dT 3' 



(13) 



just as in a disordered metal (jBeenakker and Biittikerl . 
Il992[ ). This 1/3- Fan o factor has now been confirmed ex- 
perimentally as well ( Danneau et al.l . l2007t ). 



* The T's for diffusion modes are the eigenvalues of the trans- 
mission matrix product ijt. The distribution p(T) for diffusion 
modes has a cutoff at exponential l y sma ll T ~ exp(— 2L/i), with 
/ the mean free path l|Beenakkej . Il997l ). The distribution ([T2j 
for evanescent modes has a cutoff at exp(— 47rL/A^), with A'^ 
the Fermi wave length in the heavily doped regions. In either 
case the cutoff is irrelevant for transport properties. 
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III. ANDREEV REFLECTION 



A. Electron-hole conversion 



Andreev reflection is the conversion of electron into 
hole excitations by the superconducting pair potential 
( Andreevl . Il964l ). The process is illustrated in Fig. [8] for 
the band structure of graphene. The electron excitation 
is a filled state at energy e above the Fermi energy Ep, 
while the hole excitation is an empty state at e below 
Ep- The excitation energy e is the same, so that Andreev 
reflection is an elastic process. Since electron and hole 
have opposite charge ±e, a charge of 2e is lost in the 
conversion process. This missing charge is absorbed by 
the superconductor as a Cooper pair. For e below the 
superconducting gap A electrons can enter only pairwise 
into the superconductor, and the Andreev reflected hole 
is the empty state left behind by the electron that is 
paired with the incident electron to form a Cooper pair. 

The electron and hole in Fig. [8] are taken from op- 
posite corners -iiK of the Brillouin zone, in order to 
allow the Cooper pair to carry zero total momentum. 
This corresponds to the case of s-wave pairing, common 
in conventional (low-temperature) superconductors. An- 
dreev refl^ctwniiigraphene therefore switches the val- 
leys ( Beenakkeil 120061) . The switching of valleys by 
Andreev reflection due to s-wave pairing in the super- 
conductor is analogous to the switching of spin bands 
due to singlet pairing. The latter can be detected 
by producing a spin polarizati on in the normal metal 
( De Jong and Beenakkeil . Il995l ). Analogously, the for- 
mer can be detected by producing a valley polarization 
in graphene, as we will discuss in Sec. IV. CI 

The electron and hole in Fig. [5] are both from the con- 
duction band. This intraband Andreev reflection applies 
if e < Ep. For e > Ep the hole is an empty state in 
the valence band, rather than in the conduction band. 
In undoped graphene, when Ep = 0, Andreev reflection 
is interband at all excitation energies. Interband An- 
dreev reflection does not exist in usual metals, having 
an excitation gap S> A between conduction and valence 
band. The peculiar differences between intraband and 
interband Andreev reflection are explained in the next 
subsection. 
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Andreev retro-reflection specular Andreev reflection 

FIG. 9 The left panel shows Andreev retro-reflection at the 
interface between a normal metal and a superconductor. Ar- 
rows indicate the direction of the velocity and solid or dashed 
lines distinguish whether the particle is a negatively charged 
electron (e) or a positively charged hole (ft). The right panel 
shows the specular Andreev reflection at th e interface between 
undoped graphene and a superconductor. (|Beenakkeil |2006| ) 




-kp kp 5k 
Andreev retro-reflection 



5k 

specular Andreev reflection 



FIG. 10 Dispersion relation (|14p in graphene for two values of 
the Fermi energy Ef — Tivkp, for the case of normal incidence 
[Sky — 0, Skx = Sk). Red lines indicate electron excitations 
(filled states above the Fermi level, from one valley), while 
blue lines indicate hole excitations (empty states below the 
Fermi level, from the other valley). Solid and dotted lines dis- 
tinguish the conduction and valence bands, respectively. The 
electron-hole conversion upon reflection at a superconductor 
is indicated by the arrows. Specular Andreev reflection (right 
panel) happens if an electron in the conduction band is con- 
verted into a hole in the valence band. In the usual case (left 

[ anel), electron a nd hole both lie in the conduction band. 
Beenakkeil |2006| ) 



band Andreev reflection leads to retro-reflection, while 
interband Andreev reflection leads to specular reflection. 
The linear dispersion relation ^ in graphene may be 
rewritten in terms of the excitation energy 6 — \E — Ep\^ 



e=\EF±hv{Skl + 5kiy/'^\ 



(14) 



B. Retro-reflection vs. specular reflection 



lAndreevI (1X9641) discovered that the electron-hole con- 
version at a superconductor is associated with retro- 
reflection rather than specular reflection. Retro- 
reflection means that the reflected hole retraces the path 
of the incident electron (see Fig. [9l left panel) — so 
all components of the velocity change sign. In un- 
doped graphene, in contrast, Andreev reflection is specu- 
lar (right panel) — so only th e component perp endicular 
to the interface changes sign (|Beenakkeiil2006l ). 

Inspection of the dispersion relation shows why intra- 



The ± sign distinguishes excitations in the conduction 
and in the valence band. Let the interface with the su- 
perconductor be at a; = and the electron approach the 
interface from a: > 0. Since Sky and e are conserved 
upon reflection, the reflected state is a superposition of 
the four Skx-values that solve Eq. (fH)) at given Sky and 
e. The derivative de / dSk^ is the expectation value Vx 
of the velocity in the a;-direction, so the reflected state 
contains only the two Skx-vaXnes having a positive slope. 
One of these two allowed J/caj-values is an electron ex- 
citation, the other a hole excitation. As illustrated in 
Fig. [ini the reflected hole may be either an empty state 
in the conduction band (for e < Ep) or an empty state 
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FIG. 11 Trajectories of an incident electron (red) and the 
Andreev reflected hole (blue) , for different excitation energies 
£ relative to the Fermi energy Ef, at fixed angle of incidence. 
For e < Ef the hole is in the conduction band (solid lines), 
while for e > Ef the hole is in the valence band (dashed lines) . 
The reflected trajectories rotate clockwise with increasing e, 
jumping by 180° when e — Ef- 



superconductor 




superconductor 




6 




superconductor 




superconductor 



Andreev level 



Andreev mode 



FIG. 12 The transition from retro-reflection to specular An- 
dreev reflection in a graphene channel with superconducting 
boundaries induces a transition from a localized level (left) 
to a propagating mode (right). The latter state contributes 
to thermal transport along the channel, but not to electrical 
transport. (|Titov. Ossipov. and Beenakkerll2007l ) 



in the valence band (e > Ep). A conduction-band hole 
moves opposite to its wave vector, so Vy changes sign as 
well as Vx (retro-reflection). A valence-band hole, in con- 
trast, moves in the same direction as its wave vector, so 
Vy remains unchanged and only Vx changes sign (specular 
reflection). 

The transition from retro-reflection to specular reflec- 
tion as £ increases beyond Ep is iflustrated in Fig. [TT] 
The reflection angle aout (measured relative to the nor- 
mal) first becomes greater than the angle of incidence ain, 
then jumps from -|-90° to —90° at e = Ep, and finally 
approaches — ain when e ^ Ep. 

As shown in Fig. [T^ specular Andreev reflection cre- 
ates charge-neutral propagating modes along an un- 
doped graphene channel with superc o nduct ing bound- 
aries ( Titov. Ossipov. and Beenakkeii l2007l ). In con- 
trast, retro- reflection creates bound states known as An- 
dreev levels (|Andreevl . ri964;,Kulik , .19703 . The propagat- 
ing "Andreev modes" contribute to the thermal conduc- 
tance along the graphene channel in a way which is sen- 
sitive to the superconducting phase difference across the 
channel. They may also be use d to carry a charge-neutra l 
spin current along the channel ( Greenbaum et all . l2007l ) . 
We will return to this geometry in Sec. IIII.Di when we 
consider the current across the channel (from one super- 
conductor to the other) — rather than along the channel. 



C. Dirac-Bogoliubov-de Gennes equation 

So far our discussion of Andreev reflection in graphene 
has been semiclassical, in terms of electron and hole tra- 
jectories. Quantum mechanically, the coupling of elec- 
tron and hole wave functions ^'e and ^'^ is described 
by the Dirac-Bogoliubov-de Gennes (DBdG) equation 
iBeenakker, 2006), 



fH{A) - Ep 
VA*(CTo®ro) 



A((To <E) To) 

Ep - H{-A] 



= e 



(15) 



The complex pair potential A = Age** is nonzero only 
in the superconducting region, where it couples the time- 
reversed states 



(16) 



The boundary condition for the DBdG equation at the 
edges of the graphene sheet is given by the same Eq. (fTT|) 
for both *e and ^h, 



(17) 



since we are assuming that M commutes with T. 

In the normal region A = 0, so that there ^'e and '^h 
satisfy the uncoupled equations 



(Ep- 
{Ep- 



(18) 



Andreev reflection at the normal-superconductor (NS) 
interface couples '^/^ to '^h- This coupling may be de- 
scribed by means of a boundary condition at the NS 
i nterface for the wav e function in the normal region 
(|Titov and Beenakkeii [2006 ). 



-2/3 riB cr 



(19) 



where (3 = arccos(e/Ao) G (0,7r/2) (assuming e < Aq). 
The unit vector is perpendicular to the NS interface, 
pointing from N to S. By substituting the deflnition ([T| 
of and '^fi we see that the boundary condition ^ 
couples electron excitations in one valley to hole excita- 
tions in the other valley (in accord with the description of 
Andreev reflection given in Sec. IIII.^ . In contrast, the 
boundary condition p?]) at the edges of the graphene 
sheet does not couple and '^h. 

The relation (|19p follows from the DBdG equation (fT51) 
under three assumptions characterizing an "ideal" NS in- 
terface: 

• The Fermi wave length in S is sufficiently small 
that <c; £,,Xp, where Xp = hv/Ep is the Fermi 
wave length in N and ^ = Tlv/Aq is the supercon- 
ducting coherence length. 

• The interface is smooth and impurity free on the 
scale of 
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• There is no lattice mismatch at the NS interface, so 
the honeycomb lattice of graphene is unperturbed 
at the boundary. 

The absence of lattice mismatch might be satisfied by 
depositing the superconductor on top of a heavily doped 
region of graphene. As in the c ase of a semicon- 
ductor two-dimens ional electron gas (jFagas et al 1 1200,4 



IVolkov et al.l . llQQSt ) , we expect that such an extended su- 
perconducting contact can be effectively described by a 
pair potential A in the x-y plane (even though graphene 
by itself is not superconducting). 

At normal incidence and ^'/,, are eigenstates of • 
cr, so the boundary condition implies that '^h = 
a phase factor and the electron-hole conversion happens 
with unit probability. This is entirely different from usual 
NS junctions, where Andreev reflection is suppressed at 
any angle of incidence if the Fermi wave lengths at the 
two sides of the interface are very different. 



D. Josephson junctions 

The boundary condition (|19p at a normal- 
superconducting interface depends on the phase 
of the superconductor, although this dependence is 
unobservable if there is only a single superconductor. A 
Josephson junction is a junction between two supercon- 
ductors with a controllable phase difference (/< = $i — $2- 
A current /(</>) flows from one superconductor to the 
other if ^ 0. The current flows in equilibrium, so it 
is a di ssipationles s supe rcurrent. This is the Josephson 
effect ( Josephsonl Il964l ). Since / is 27r-periodic in (f>, 
there exists a maximal supercurrent Ic that can flow 
between the superconductors. This is called the critical 
current of the Josephson junction. 

There is a thermodynamic relation ( AndersonL 1 1 963^ 




FIG. 13 Josephson junction, formed by a graphene layer (N) 
with two superconducting electrodes ( S) a distance L apart , 
havin g a phase difference </) = $i — $2. (|Titov and Beenakkerl 
I2OO6I ) 
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FIG. 14 Critical current Ic and IcRn product of a bal- 
listic Josephson junction (length L short compared to the 
width W and superconducting coherence length as a func- 
tion of the Fermi energy Ef in the normal region. Small 
and large \Ef\ asymptotes are indicated by dashed lines. 



(20) (|Titov and Beenakke: 



optotes J 

Ellooi) 



between the supercurrent / and the derivative of the free 
energy F with respect to the superconducting phase dif- 
ference. The free energy can in turn be related to the 
excitation spectrum, which itself follows from the DBdG 
equation. At zero temperature and in the short-junction 
limit (separation L of the two NS interfaces <C £) th e 
resulting relation is (Beenakker and Van Houtenl . Il992() 



(21) 



with En < Ao the energy of a (spin-degenerate) bound 
state in the Josephson junction. 

To calculate the supercurrent (in zero magnetic field) 
one therefore needs to solve the two eigenvalue equa- 
tions ([TS]) (with A = 0) in the strip < x < L, 
<y <W (see Fig. [T3]). At x = 0,L there is the phase- 
dependent boundary condition p9p , which couples to 



^'h, while the boundary condition (fTT|) at y = 0, is 

phase-independent and does not c ouple to '^h- 

Th e result of this calculation (jTitov and Beenakkeii 
I2OO6D is that the critical current is given, up to numerical 
coefficients of order unity, by 



eAo 



max(H//L, VF/Aj.). 



(22) 



(The dependence on the boundary condition at y = 0, VF 
can be neglected under the assumption L PF of a 
short and wide junction.) At the Dirac point Ep ~ 
one has Aj? — > cxd, so the critical current reaches its min- 
imal value ~ (eAo/h) x W/L (see Fig.[T31 upper panel). 
Instead of being independent of the length L of the junc- 
tion, as expected for a short ballistic Josephson junction, 
the critical current at the Dirac point has the diffusion- 
like scaling cx 1/L. This is another manifestation of the 
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FIG. 15 Product of the critical current 7c and the normal 
state resistance Rn versus gate voltage V^ate, measured at 
T = 30 mK in the Josephson junction of Fig. [2] The carrier 
density in the graphene layer is linearly proportional to Vgato, 
while the Fermi energy Ef (x v'V'gato. The resistance Rm is 
measured in the presence of a sma ll magnetic field to d rive 
the electrodes in the normal state. (jHeersche et al.l . l2007l ) 



"pseudo-diffusive" dynamics discussed in Se c. III. Dl 

Since the nor mal-state resistance scales as ( Katsnelsonl . 
l2006bl : iTwo^dlo et al.. . ,200(1 



l/i?Ar = G 



N 



max{W/L,W/XF), 



(23) 



the theory predicts that the IqRn product remains of 
order Ag/e (up to a numerical prefactor) as the Fermi 
level passes through the Dirac point (Fig. [HI lower 
panel). The experimental result of lHeersche et al.1 ( 20071 ) 
for the Josephson junction of Fig. [2l shown in Fig. [TSl 
is qualitatively similar to the theoretical prediction, but 
there are significant quantitative differences: The ex- 
perimental IcRn product at the Dirac point is about 
60 /i-V w 0.5Ao/e, more than twice the theoretical pre- 
diction, and the increase at higher carrier densities is 
much larger than predicted. It is quite likely that disor- 
der in the experimental sample, which is not included 
i n the calculation, is responsib le for these differences 
(|Du. Skachko. and Andreil . 120071) . 



E. Further reading 

In the spirit of a Colloquium, we have only discussed 
the basic physics of Andreev reflection in graphene. In 
this subsection we give some pointers to the literature on 
other aspects of this topic. 

The pseudo-diffusive dynamics at the Dirac point, dis- 
cussed in Sec. lIII.Dl in connection with the critical current 
Ic of an SNS junction, extends to t he entire current-phase 
relati onship I{(f>) in equilibrium ( Titov and Beenakkeii 
I2006D . a s well as to the dissipatiye curr ent out of equi- 
librium ( Cuevas and Lew Yevatil . l2006l ). In each case. 



a short and wide strip of undoped ballistic graphene 
(length L short compared to width W and superconduct- 
ing coherence length ^)^ behaves as a disordered metal 
having the same normal-state conductance Gn ■ 

Pseudo-diffusive dynamics also governs the con- 
ductance Gns through a ballistic graphene strip 
(L «C W) having a sin g le sup e rconducting contact 
lAkhmerov and Beenakkol l2007at iPrada et al.l [2007), 
in the sense that the ratio Gns/Gn at the Dirac point 
is the same as for a disordered metal. The correspon- 
dence holds only for voltages V small compared to hv/L. 
At larger voltages the current-voltage characteristic of a 
ballistic NS juncti o n in graphene has unusual f e atures 
(iBeenakkeil 12006 : Bhattachariee arid Senguptal . 120061 : 
iBhattachariee. Maiti. and Sengupta . 20071 ) — "without a 
diffusiv e analogue. These h ave been studied expe rimen- 
tally bv lShailos et~an (|2007t) and lMiao etld] (|2007f ). Sim- 
ilarly unusual I- V characteristic s have been predicted in 
bilayer graphene ( Ludwid . l2007t ). 



The Dirac-Bogoliubov-de Gennes equation of Sec. lIII.Cl 
assumes isotropic (s-wave) pairing in the superconduc- 
tor. The equation may be readily modified to the case of 
anisotropic (c?-wave) pairing, relevant for NS contacts be- 
tween graphene and a high-temperature superconductor. 
The conductance in the two cases has been compared by 
iLinder an d Sudb^l ([2007') . Another modification, studied 
by Wchlin g et al.l ([2007), is to include electrical or mag- 
netic scattering potentials in the superconducting region. 

More exot ic {px + iPv or d^,2 ,,2 + idxv) pairings may 
be p ossible ( Jiang et al.l . 120071 : lUchoa and Castro Netol . 
|2007| ) if graphene could be chemically modified to be- 
come an intrinsic superconductor (rather than having 
the superconductivity indu ced by the proximity effect). 
IChaemi and Wilczekl (|2007l ) have argued that the special 
topological properties (nonabelian statistics) of vortices 
in a superconductor with px -f ipy pairing apply as well to 
the s-wave DBdG equation psp , if the superconductivity 
can be induced in undoped graphene. 

The idealized model of the NS interface discussed in 
Sec. nil. CI can be much improved, in particular to in- 
clude the effects of lattice mismatch and a selfconsistent 
calc ulation of the induc ed pair potential . Some iiumeri - 
cal (|Wakabavashil l2003f ) and analytical (|Tkachovl |2007| ) 
work goes in this direction. 



^ This short-junction Umit is essential: Pseudo-diffusive dynam- 
ics in SNS junctions breaks down if L becomes larger t h an W 
l lGonzalez and Perfettol . |2007| : iMoehaddam and Zarevanl. 120061) 
or if L becomes larger than e llTitov. Ossipov. and Beenakken . 
I2OO7I) . A tunnel barrier l|Ma iti anT'Seng uptan2007t) or r>-n iunc- 
tion l|OssiDov. Titov. and B ccnakkcr , 2003) in the normal region 
also spoils the pseudo-diffusive analogy. 
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FIG. 16 Classical trajectories of an electron in the presence 
of a uniform electric field in the a;-direction. All three trajec- 
tories are at the same energy, only the component Py of the 
momentum transverse to the field lines is varied. The two 
black trajectories are for py > 0, while the green trajectory is 
for Py =0. The electron is in the conduction band of graphene 
for X < Q (solid trajectories, velocity parallel to momentum) 
and in the valence band for j: > (dashed trajectories, veloc- 
ity antiparallel to momentum). Solid and dashed trajectories 
are coupled by Klein tunneling. 




IV. KLEIN TUNNELING 




A. Absence of backscattering 

The massless carriers in graphene respond quite differ- 
ently to an electric field than ordinary massive carriers. 
Because the magnitude v of the carrier velocity is inde- 
pendent of the energy, an electron moving along the field 
lines cannot be backscattered — since that would require 
w = at the turning point. The absence of backscatter - 
ing was discovered bv lAndo. Nakanishi. and Saitol ( 19981 ) 
in carbon nanotubes, where it is responsible for the 
high conductivity in the presence of disorder. The two- 
dimensional nature of the dynamics in graphene intro- 
duces some new aspects. 

Electron trajectories in the linear electrostatic poten- 
tial U {x) = Tx are shown in Fig. [111 The trajectories are 
deflected by the electric field for py ^ 0, but for py = Q 
no backscattering occurs. The electron is able to propa- 
gate through an infinitely high potential barrier because 
it makes a transition from the conduction band to the 
valence band (see Fig. I17p . In this transition its dynam- 
ics changes from electron-like to hole-like in the following 
sense: 

The equation of motion 



dr 
HI 



dE_ 
dp 



9 

V p 

E-U' 



(24) 



at energy E with — [E — C/)^, implies that the 

velocity v = dr/dt of the electron is parallel to the mo- 
mentum when it is in the conduction band (U < E) 
and antiparallel when it is in the valence band (U > E). 
States with v parallel to p are called electron-like and 
states with v antiparallel to p are called hole-like. By 
making the transition from electron-like to hole-like dy- 
namics, the electron can continue to move in the same 
direction even as its momentum along the field lines goes 
through zero and changes sign. 

In classical mechanics, backscattering is only avoided 
for Py = (so only if the electron moves along 



FIG. 17 Band structure of a single valley at two sides of a 
potential step (height Uq, width d). The equilibrium Fermi 
energy Ef is the same at both sides, so that for Uo > Ef an 
electron just above the Fermi level is in the conduction band 
at one side and in the valence band at the other side. Blue 
arrows indicate the electron velocity, which is parallel to the 
wave vector (or momentum) in the conduction band (left) and 
antiparallel in the valence band (right). 



the field lines). In quantum mechanics an elec- 
tron can tunnel from the conduction into the va- 
lence band, thereby avoiding backscattering, also for 
a small but nonzero Py. Such tunneling from an 
electron-lik e to a hole-like s tate i s called interband 
tunn e ling (lAronov and Pikusl. 1967t Kane and Blountl . 



1969|i^ IWeiler. Zawadzki. and Laxl.ll967l) or K lein tunnel- 
ing (Katsn elson. Novoselov. and Geiml . l2006l ). because of 
an analogou s effect in relativistic quantum mechanics 
(iKleinl . 119291) . 1° 



The probability of Klein tunneling of a relativistic 
electro n in a uniform electric field was calculated by 
ISauteil (Il93ll) — with an exponentially small result 
due to the finite electron mass. The case of mass- 
less particles, relevant for graphene, was considered by 
ICheianov and Fal'kol (|2006D . Pairs of electron-like and 
hole-like trajectories at the same E and Py have turn- 
ing points at a distance dmin = 2w|py|/jF. The tunneling 
probability has an exponential dependence on this sepa- 



Klein t unneling is considered paradox ical in the relativistic 
context llCaloeeracos and Dombevl . Il999t) . because the hole-like 
states into which the electron tunnels are unphysical antiparticle 
states in the Dirac sea. There is no paradox in the context of 
graphene, where the role of the Dirac sea is played by the valence 
band (see Fig. (Tjl. 
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FIG. 18 n-p-n junction in graphene: a) cross-sectional view 
of the device, b) electrostatic potential profile U{x) along 
the cross-section of a). The combination of a positive volt- 
age on the back gate and a negative voltage on the top gate 
produces a central p-doped region flanked by two n-doped 
regions, c) Optical image of the device. The barely visible 
graphene flake is outlined with a dashed line a nd the dielec- 
tric l ayer of PMMA appears as a blue shadow. (jHuard et al.l . 
I2007l 1 



ration, 

T{py) = exp(-7r|p.y|d,nin/2?i) = eK.p{~TTvpy/hT), (25) 

provided that the longitudinal momentum p^^ at a; — > 
— oo and at a; ^ oo (where the electric field is as- 
sumed to vanish) is sufhciently large: 



\p%\pT\^\PylVf^- 



(26) 



Transmission resonances occur when a p-n interface 
is in series with an n-p interface, forming a p-n-p or 
n-p-n junctio n fe atsnelson. Novoselov. and Geim . 



2006 



Milton Pereira et al 



Milton Pereira. Vasilopoulos. and Peetersl 
Silvestrov and Efetovl . l2007t) . The unit 
sion at py 



2006; 



2007 



transmis- 

forbids transmission resonances at normal 



incidence, in marked contrast with conventional resonant 
tunneling through a double-barrier junction. 



B. Bipolar junctions 

Klein tunneling is the mechanism for electrical conduc- 
tion through the interface between p-doped and n-doped 
graphene. Su c h a bi polar junction is illustrated in Fig.[T51 
( Huard et al.1 . |2007[ ). A top gate creates an electrostatic 
potential barrier, so that the Fermi level lies in the va- 
lence band inside the barrier (p-doped region) and in the 
conduction band outside the barrier (n-doped region). 
The carrier density n-carrior is the same in the n and p re- 
gions when the Fermi energy Ep is half the barrier height 
Uq. In this case the Fermi momenta pp = ^kp in both n 
and p regions are given by pp = Uo/2v = J-d/2v, with d 
the width of the n-p and p-n interfaces and = Uo/d the 
electric field (up to a factor of electron charge) in that 
interface region. 

The width d is of the order of the separation b e tween 
the graphene layer and the top gate. iHuard et al. I (|2007D 
estimate d « 80 nm for their device. The Fermi wave vec- 



tor kp — y/ji 



> 10-- 

' cm" 



, for typical carrier 

densities OI ncarrier <L ""^^ 

ditions kpd > 1 the p-n and n-p interfaces are smooth 
on the scale of the Fermi wave length. This is the regime 
of applicability of the expression (pS)) for the Klein tun- 
neling probability, since the condition of large lon- 
gitudinal momentum can be rewritten as 



Pp > \py\,h/d. 



(27) 



The conductance Gp_„ of a p-n interface follows by 



integration of Eg. (|25D over the transverse momenta, with 
the result ([Ch eianov and Farkg . ,2006. ) 



4e^ W 



dPyT{Py)^—-xl-, (28) 



The factor of 4 accounts for the twofold spin and valley 
degeneracy and W is the transverse dimension of the in- 
terface. The integration range may be extended to ±00 
because T(py) is vanishingly small for \py\ larger than 
^JhTjv ~ ppj \Jkpd pp. Notice that Gp-„ is smaller 
than the baUistic conductance GbaiUstic — {^e^ /h)kpW/T: 
by the same factor y/kpd that characterizes the smooth- 
ness of the interface. 



The asym ptotic res ult ||25| I. derived by ICheianov and FaT 
ll2006h and lAndreevI {2007), shoul d foll ow from the general 
Kummer- function formula of .Sauterl lll93lh upon substitution of 
the electron mass m by Py/v and taking the limit H26p . The 
asymptotic limit taken by Sauter corresponds to the opposite 
regime \py\ 3> \Px^\, |pS"*|> yJ'hT jv in which T(py) is exponen- 
tially small. 



This assumption of a constant electric field in the interface re- 
gion requires perfect screening by the carriers in graphene of the 
electric field produced by the gate. The lack of screening at the 
p-n interface due to the vanishing carrier density enhances the 
local electric field by a factor {e^kpd/ kHv)^/^ , with k the dielec- 
tric constant llZhang and FogleH. 1200% . The value of ft can be as 
low as 2.4 for a Si02 substrate and as large as 80 for graphene 
on water. 
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FIG. 19 Two trajectories in a p-n-p junction, the lower one 
(transmitted) in zero magnetic field and the upper one (re- 
flected) in a small but nonzero field. Because only trajectories 
with an angle of order 1 / \/kFd, <^ 1 around normal incidence 
are transmitted through the p-n and n-p interfaces, a rela- 
tively weak magnetic field suppresses the series conductance 
of the inter faces by bending the trajecto ries away from normal 
incidence. (|Cheianov and Fal'kol . [20061 ) 



p and n regions propagate parallel t o th e in- 
terface (iLukose. Shankar. and Baskara R 120071: 

iMilton Pereira. Peeters. and Vasilopoulosl 20071 ). The 
corresponding classical trajectories are illustrated 
in Fig. [20l The direction of propagation along the 
interface is the same in both the p and n regions 
( Abanin and Levitovl . |200"7[) — while the direction of ro- 
tation is opposite. The snake-shaped trajectory centered 
at the interface has a mixed electron-hole character. 
This is the ambipolar analogue of the snake states that 
are known to exist in a nonhomogeneous magn e tic fie ld 
( Ghosh et al. I, l2007t iMiilled . Il992l : iRakvta et all . I2007D . 

The conductance in the high-field regime B > B^, is not 
fully suppressed, but it no longer scales with the width 
W of the junction. We will calculate it in Sec. IV. CI 



D. Further reading 



y 

q'^Vj' ' 

FIG. 20 Electron trajectories along a p-n interface in a mag- 
netic field B > B, (when there is no transmission through the 
interface). The electron rotates in opposite directions in the 
conduction band (solid trajectories) and in the valence band 
(dashed). The trajectory centered at the interface (green) 
represents an "ambipolar snake state". 



C. Magnetic field effects 



[Ch eianov and Fal'kol (|2006l ) have predicted that a rela- 
tively weak magnetic field B ~ (h / e) y/kp/Ti^ will sup- 
press the conductance of an n-p-n or p-n-p junction (of 
length L) below the series conductance of the individual 
interfaces, as a consequence of the strong angular depen- 
dence of the transmission probability (|25p . The mecha- 
nism is illustra ted in Fig. [TOl The effect is not observed 
in the device o f lHuard et al .' ('2007') , presumably because 
of disorder. (See iFogler et al... (2007,) for a calculation of 
the conditions required for ballistic transport, which are 
only marginally met in existing experiments.) 

The conductance of a single p-n interface becomes 
magnetic field dependent on the much larger field scale 
B^ = T I ev ~ {fi/ e)kp / d at which the cyclotron radius 
^cyci = hkp/E B becomes comparable t o the w idth d of 
the interface. IShvtov. Gu. and" Levitovl (|2007l) have cal- 
culated that the angle of incidence ^max which is maxi- 
mally transmitted rotates away from normal incidence to 
a value ^max = ± arcsin(i3/i3*). The effect on the con- 
ductance (P5)l of the p-n interface is a suppression by a 
factor [1- (B/S*)2]3/4. 

For _B > no transmission is possible through 
the p-n interface. Instead, states in both the 



Just as in Sec. IIII.EI we mention here some papers for 
further reading on this topic. 

Klein tunneling in a carbon bilayer differs fun- 

damentally from Klein tunnelin g in a monolayer 

([Katsnelson . Novos elov. and Geiml . I2006D . The bilayer 



still has a gapless spectrum [i n the absence of a po t ential 
difference between the layers ( McCann and Fal'kol . l2006t 
lOhta et "all |2006[ )]. so interband tunneling can happen 
with high probability. However, at normal incidence the 
probability is — while it is 1 in the monolayer. Al- 
though electrons in a bilayer are not massless, as they 
are in a monolayer, they still have a definite chirality (di- 
rection of motion tied to direction of pseudospin). Klein 
tunneling in a carbon bilayer is therefore different from 
interband tunneling in a gapless semiconductor. For ex- 
ample, the chirality forbids transmission resonances at 
normal incidence. 

The perfect transmission at normal incidence in a 
monolayer is a robust effect with regard to the shape of 
the electrostatic potential profile at the p-n interface (all 
that is needed is a potential which is smooth on the scale 
of the lattice constant). A time dependent electric field 
parallel to the interface, however, can suppress the trans- 
mission even at normal incidence (iFistul and Efetovl . 



l2007t iTrauzettel. Blanter. and Morpurgol l2007t ). The 
suppression is strongest if the frequency u of the radi- 
ation satisfies the resonance condition uj = 2v\p\/h at 
some point in the interface region. 

Bipolar junctions may appear naturally in disordered 
graphene, when the random electrostatic potential land- 
scape produ ces alternating regi ons of p-type and n- 



type doping ( Martin et al.l . l2007l ). Classical percolation 



through such a r andom network o f bipo lar junctions has 
been studied by ICheianov et al. 1 IIOO^)- At zero Fermi 
energy (when the areas of p-type and n-type doping are 
the same), the percolation length remains finite because 
of Klein tunneling. 

Klein tunneling is also responsible for the finite life 
time of an electron state bound to a charged impurity in 
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graphene ( Shvtov. Katsnelson. and Levitovl . [20071 ) . Such 
quasi-bound states exist for /? = Ze^/nfiVF > 1/2, with 
Ze the impurity charge. The discrete states exist in the 
conduction band near the impurity, but they are cou- 
pled by Klein tunneling to the continuum of states in 
the valence band away from the impurity. The resulting 
resonances (having width-to-energy ratio e~'^^^) may be 
observable by measuring the local density of states with 
a scanning probe. 



V. ANALOGIES 

In the previous two sections we have discussed NS and 
p-n junctions separately. In this section we address the 
analogies between these two systems. Both involve the 
coupling of electron-like and hole-like states, either by 
the superconducting pair potential (in the NS junction) 
or by the electrostatic potential (in the p-n junction). 
An obvious difference is that the two types of states 
lie at the same side of the NS interface but at oppo- 
site sides of the p-n interface. The analogies, therefore, 
i nvolve a reflection of t he geometry along the interface 
(iTworzvdlo et al.l . l2007t) . 



A. Mapping between NS and p-n junction 



A precise mapping (jBeenakker et al.l . l2007l ) between 
NS and p-n junctions is possible under two conditions: 

• The electrostatic potential U in the p-n junction is 
antisymmetric, U{—x,y) — —U{x,y), with respect 
to the p-n interface at a; = 0. 

• The NS interface may be described by the boundary 
condition (fT9t at a; = 0. 



A uniform perpendicular magnetic field B may or may 
not be present. Under these conditions a p-n junction 
has the same excitation spectrum as an NS junction for 
Ep — and excitation energies e ^ Aq. 

This correspondence follows from the fact that, if 
is an eigenstate of the Dirac Hamiltonian ^ of the p-n 
junction with eigenvalue^^ e, then we can construct an 
eigenstate (^'e, 5'/;) in the normal part a; > of the NS 
junction by 



*e(a;,y) = *(a;,j/), 

*;,(a;, y) = ie'^^ia^ ® To)*(-a;, y) = 7'*(x, y). 



(29) 



Here V = ie^**(cra; <8i tq)TI with TZ the reflection oper- 
ator {x 1-^ —x). Since 5* is continuous at x = 0, the 
boundary condition at the NS interface is automat- 
ically satisfied for e <C Aq. Furthermore, from i/vj/ = evj/ 



and VH{A) = -H{-A)V (with A = Bxy) it follows 
that ^'e and '^h satisfy the DBdG equation in the 
normal region. 

The applicability of the mapping may be extended in 
several ways: The p-n junction may have boundaries de- 
scribed by the boundary condition 5'(r) = M(r)^'(r) for 
r at the boundary. (We may assume that this relation 
holds for all r by setting M = 1 when r is not at the 
boundary.) The mapping to an NS junction still holds, 
provided that M commutes with V , which requires 



(30) 



{<Jx ® To)M{x, y) = M{-x, y){ax ® tq). 



This ensures that the transformed wave function ((29|) 
in the NS jimction satisfies the corresponding bound- 
ary condition \n\ . For example, an armchair boundary 
along the x-axis (with M cx (j^ independent of x) satisfies 
the requirement pO|) , but a zigzag boundary along the x- 
axis (M oc az) does not. A pair of zigzag boundaries at 
x = ±L (with M{±L, y) — zbcr^ (8)i/z), on the other hand, 
do satisfy the requirement ([50]) . 

The Dirac Hamiltonian ^ of the p-n junction may 
also contain an additional term SH without spoiling the 
mapping to the NS junction, provided that 6H anticom- 
mutes with the operator product VT, 



VTSH = -6HVT. 



(31) 



Since the spectrum of the p-n junction is symmetric around zero 
energy, it suffices to consider energies £ > 0. 



Considering the two examples of a SH mentioned in 
Sec. III.Bl we see that the mass term preserves the map- 
ping if fi{x,y) — fi(—x,y), while the valley-dependent 
vector potential should satisfy Ax{—x,y) = —Ax(x,y), 



B. Retro-reflection vs. negative refraction 

We apply the mapping of the previous subsection to an 
abrupt p-n junction, as shown in Fig. 1211 With "abrupt" 
we mean that the width d of the potential step at the 
p-n interface should be small compared to the Fermi 
wave length Xp = hv/Uo- In Sec. IIV.BI we discussed 
the opposite regime d^ Xp of a smooth interface, when 
only electrons approaching the interface near normal in- 
cidence are transmitted. For an abrupt interface the 
transmission probability is large also away from normal 
incidence j and an tuiusual effect of negative re fraction 
appears ( Cheianov. Fal'ko. and Altshuled . HOO^): Upon 
crossing the p-n interface the sign of the tangential ve- 
locity component is inverted. 

The lower panels in Fig. [^T] show how the classical 
trajectories in the p-n and NS junctions are mapped 
onto each other by reflection along the interface. Retro- 
reflection in the NS junction (inversion of the tangen- 
tial velocity component upon conversion from electron 
to hole) maps onto negative refraction in the p-n junc- 
tion. As the excitation energy e increases beyond the 
step height Uq, negative refraction crosses over into pos- 
itive refraction at the p-n junction in the same way that 
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FIG. 21 Comparison of two systems that can be mapped onto 
each other by the transformation (|29p . The upper graphs 
show the electrostatic potential profile (solid lines) of a p-n 
junction (left) and of the corresponding NS junction (right, 
with Uao C^o)- The upper right graph also shows the su- 
perconducting pair potential A (dashed line). The excitation 
spectrum of the two systems is the same for e <^ Aq. Clas- 
sical trajectories in the two systems are related by reflection 
along the interface, as shown in the lower graphs for e = 
(solid lines indicate the electron-like trajectories and dashed 
lines the hole-like trajectories). 



retro-rcflcction crosses over into specular reflection at the 
NS junction (compare Figs. [TT] and [^^ . 

Because the mapping (I^Hl) is quantum mechanical, not 
only the trajectories are mapped onto each other but the 
full diffraction pattern together with the quantum me- 
chanical transmission and reflection probabilities. For 
example, when e ^ Uq, Ag the NS junction has a prob- 
ability Reh — cos^ 6 for Andreev reflection (electron 
to hole) and a probability R^e — sin^ 6 for normal re- 
flecti on (electron to electron). This result (jBeenakkeii 
hood ) agrees with the transmission and reflection prob- 
abilities T,R in an a brupt p-n junction calculated by 
[Ch eianov and Fal 13 (|2006D . upon mapping Reh ^ T 
and i?ee 1-^ R- For normal incidence {9 = 0) both retro- 
reflection and negative refraction happen with unit prob- 
ability. 

Negative refraction was first discovered in optics,^'* 
where it is used as a way to make a flat lens known 
as a Ve selago lens jVescl a gQ, 1968() . [For a tuto- 
rial, see iPendrv and Smithl ( 20041 ).] As calculated by 



^'^ The most direct analogy is with the work of iNotomil ||2000| '| on 
negative refraction in two-dimensional photonic crystals with the 
same honeycomb lattice as graphene. 




e<Uo 



FIG. 22 Trajectories of an incident and refracted electron at a 
p-n interface, for different excitation energies e relative to the 
potential step height Uo, at fixed angle of incidence. For e < 
Uo the refracted electron is in the valence band (dashed lines), 
while for e > Uo it is in the conduction band (solid lines). The 
refracted trajectories rotate counter-clockwise with increasing 
£, jumping by 180° when e — Uo- The transformation x —x 
maps this transition from negative to positive refraction onto 
the transition from retro-reflection to specular reflection in 
the NS junction of Fig. \TT\ 



n 



P 



n 



FIG. 23 Classical trajectories (dotted lines) in an n-p-n junc- 
tion at an energy e = that is halfway the potential step 
across the n-p and p-n interfaces, so that the refraction pre- 
cisely inverts the angle of incidence. A scatterer in the n- 
region (blue horizontal line) has an inverted image in the cen- 
tral p-region and a noninverted image in the other n-region. 
This is the principle of operation of the Veselago lens. 



[Ch eianov. Fal'ko. and Altshuleil (|2007l ). an abrupt p-n 
interface produces an inverted image in the p-region of 
a scatterer in the n-region. An n-p-n or p-n-p junc- 
tion inverts the image twice, reproducing the original 
image at the other side of the junction (see Fig. [^5)1 . 
The Veselago lens in graphene is not ideal: Negative re- 
fraction only produces a perfect focus at £ = 0, while 
at other energies the focus is spread into a caustic. 
Caustics (focal lines, rather than focal points) also ap- 
pear if the p-n interface is curved rather than straight 
(|Cserti. Palvi. and Peterfalvill2007t) . 
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FIG. 24 Schematic top view of a graphene nanoribbon con- 
taining an interface between a p-doped and n-doped region 
(left panel) and between a normal (N) and superconduct- 
ing (S) region (right panel). Electron- like and hole-like edge 
states in the lowest Landau level are indicated by solid and 
dashed lines, respecti vely, with arrows pointing in the direc- 
tion of propagation. (|Tworzvdlo et af] , 120071 ) 



C. Valley-isospin dependent quantum Hall effect 



In Sec. III. CI we mentioned that the edge states in the 
lowest Landau level are valley polarized, with a valley 
isopin u determined by the boundary condition pT|) at 
the edge. Here we discuss how this valley polarization 
can be measured in a conduction experiment on either a 
p-n junction or an NS junction. 

The two geometries are compared in Fig. [531 Electron- 
like and hole-like valley-polarized edge states hybridize 
along the p-n or NS interface to form a valley-degenerate 
electron-hole state. (In the p-n case, this state corre- 
sponds classically to the snake-shaped trajectory in Fig. 
[501) The two-terminal conductance G = GoTeh is deter- 
mined by the probability Teh that an electron-like state is 
converted into a hole-like state at the opposite edge (with 
Go = 2e^ /h in the p-n juncti on and Gp = 4e^//t in the 
NS jun ction) . As shown by lAkhmerov and Beenakkerj 
(l2007bD and iTworzvdlo et all (12001 7 in the absence of 
intervalley scattering this probability 



Teh = i(l -COS$) 



(32) 



depends only on the angle $ between the valley isospins 
of the electron-like state at the two edges. 

Eq. (|551) assumes that the electron-like and hole-like 
edge channels at one edge have opposite valley isospins 
(ii'i for the left edge and ±i//j for the right edge).^^ 
Since the unidirectional motion of the edge states pre- 
vents reflections, the total transmission matrix itotai = 



One factor of two in Go comes from the spin degeneracy. The NS 
junction has one more factor of two because the electron-to-hole 
conversion transfers two electrons across the junction. 

This is generally the case, with one exception; A p-n junction in 
a zigzag nanoribbon has electr on-like and hole- l i ke edg e channels 
with identical valley isospins iTworzvdlo et ahl l2007t) . 
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FIG. 25 Conductance of an armchair nanoribbon containing 
the potential step U[x) — i [tanh(2a;/Z/) + l]Uac, calculated 
numerically from a tight-binding model in a perpendicular 
magnetic field (magnetic length Im = ^Ti/eB = 5 a). The 
step height C/oo is varied from below Ef (unipolar regime) 
to above Ef (bipolar regime), at fixed Ef = Tiv/lm and 
1/ = 50 a. The solid curves are without disorder, while the 
dashed curves are for a random electrostatic potential land- 
scape (correlation length ^ = 10 a). The different colors corre- 
spond to a different number N of hexagons across the ribbon, 
and hence a different width W ^ {M + 3/2)a: A/" = 97 (red 
curves), 98 (blue), and 99 (green). T he dashed horizontal line 
marks the plateau at G = i x 2e^//i. l|Tworzvdlo et al.l ■ [2007^ 



tRtinttL from one edge to the other edge is the product 
of three 2x2 unitary matrices: the transmission matrix 
from the left edge to the p-n or NS interface, the trans- 
mission matrix tint along the interface, and the transmis- 
sion matrix tn from the interface to the right edge. In the 
absence of intervalley scattering tint = e*'^'"'ro is propor- 
tional to the unit matrix in the valley degree of freedom, 
while 



tx 



J4>> 



'^x){+i^x\+e'^-\-ux){-iyx\ (33) 



(with X ~ L,R) is diagonal in the basis | ± i^x) of eigen- 



states of Vx t- The phase shifts 



need not be 



determined. Evaluation of the transmission probability 

Tek ^ \{ + l^L\ttot.l\ ~ l^R)f (34) 

leads to the conductance 

G= iGo(l-cos$), (35) 

with cos<I> = vl ■ i^R- 

The angle $ = AttW / 3a -\- tt between the valley isospins 
at two opposite armchair edges depends on the width W 
(as defined in Fig. [3]): $ = vr if 2W/a is a multiple of 3, 
|$| = 7r/3 if it is not (see Fig. [5]). A tight- binding model 
calculation of an armchair nanoribbon containing a po- 
tential step (Fig. [55|l indeed shows that the conductance 
as a function of the step height switches from a plateau 
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FIG. 26 Experimental conductance of a gate-controlled p-n 
junction in graphene. The conductance of the n-doped re- 
gion at one side of the interface is fixed at \fi\e^/h, with 
/i = 2, while the conductance \f2\e^/h at the other side 
of the interface is varied by the gate voltage (values of /2 
are indicated, with negative numbers corresponding to a p- 
doped region). In the unipolar regime (/1/2 > 0) the con- 
ductance of the junction is given by G = min(|/i|, |/2|)e^/ft, 
while in the bipolar regime (/1/2 < 0) the conductance is 
the Ohmic series conductance G x h/e ^ = l/iM/d/il + IM). 
(jWilliams. DiCarlo. and Marcud [2OO7I ) 



at the ^-independent Hall conductance Gq in the unipo- 
lar regime {n-n junction) to a ^-dependent value given 
by Eq. (j35|) in the bipolar regime {p-n junction). The 
plateau persists in the presence of a random potential, 
provided it is smooth on the scale of the lattice constant 
so no intervalley scattering is introduced. 

The valley-isospin dependence of the quantum Hall ef- 
fect makes it possible to use strain as a means of variation 
of the height of the conductance plateaus. As mentioned 
in Sec. III.Bl strain introduces a valley-dependent vec- 
tor potential in the Dirac equation, corresponding to a 
fictitious magnetic field of opposite sign in the two val- 
leys. This field rotates the Bloch vector of the valley 
isospin around the 2;-axis, which in the case of an arm- 
chair nanoribbon corresponds to a rotation of the valley 
isospin in the x-y plane. 

In the high-magnetic field exper iments of 
IWilliams. DiCarlo. and Marcud (120071 ) and 

lOzvilmaz et~al1 (|2003) the p-n junction has a quantized 



con ductance, see Fig 



by lAbanin and Levito\r^2007f ) 

conductance Gsorics = GpGn/ (Gp 



This 



has been explained 
1 the Ohmic series 
Gn) of the quantum 



Hall conductances Gp,Gn in the p-dopcd and n-doped 
regions (each an odd multiple of the conductance 
quantum 2e'^/h). Ohm's law would apply if the system 
is sufficiently large that a local equilibrium is established 
at the interface, while the non-Ohmic result would 
be expected for smaller systems. 
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FIG. 27 Plot of the density of states p(e) for the p-n junction 
shown in (a). The dotted line is the value in the isolated 
p and n regions, which is energy independent for |e| <^ Uo- 
The density of states vanishes at the Fermi level (e — 0), 
according to Eq. (|36p . The NS junction shown in (b) has the 
same density of states. In both the NS and p-n geometries 
the suppression of the density of states is due to destructive 
interference of the electron-like and hole-like segments of the 
periodic orbits at the Fermi level (indicated by solid blue and 
dashed green trajectories). jBeenakker et all . 120071 ) 



D. Pseudo-superconductivity 

The correspondence between NS and p-n junctions of 
Sec. IV. Al implies that quantum effects associated with 
superconductivity, such as the proximity effect and the 
Josephson effect, have analogues in nons uperconducting 
bipolar graphene ( Beenakker et al.ll2007( ). 

Such "pseudo-superconductivity" is demonstrated in 
Fig. [571 which plots the density of states p{e) in a p-n 
junction with an abrupt interface. The p and n regions 
have the same Fermi energy Uq and zigzag boundaries at 
X = ztW. The width W is assumed to be large compared 
to the Fermi wave length Xp — hv/Uo- The density of 
states, smoothed over rapid oscillations, vanishes linearly 
as 



p{e) = po\e\/ET 



(36) 



for small with po = (2[/o/7r)(?iu)~^ the density of 
states (per unit area and including spin plus valley de- 
generacies) in the separate p and n regions. The en- 
ergy Et = fiv/2W is the Thouless energy (which is 
^ p, for W ^ Xp). This suppression of the density 
of states at the Fermi level by a factor e/Ex is pre- 
cisely analogous to an NS junction, where the density 
of states is suppressed by the superconducting proximity 
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FIG. 28 Persistent current through a ring containing an 
abrupt p-n interf ace, as a function of the magnetic flux 
through the ring. (jBeenakker et all , |2007| ) 



listic Josephson junction. 

This concludes our review of Andreev reflection 
and Klein tunneling in graphene. The analogies dis- 
cussed in this Section will hopefully be validated in 
the near future by ongoing experiments on bipolar 
junctions and Josephson junctions. From a differ- 
ent perspective, the correspondence between these 
two phenomena offers the intriguing opportunity to 
observe superconducting analogies in non-electronic 
systems governed by the same Dirac equation as 
graphene. An example would be a two-dimensional 
photonic crystal on a honeycomb or triangular lat- 
tice (iGarcia-Pomar. Cor tiio. a nd Nieto-Vesperinas , 



2007; 



Raghu and I faldam 



Sepkhanov. Bazaliv. and Beenakkerl 1200 



2006 : 



effect ( Titov. Ossipov. and Beenakkeij . l2007[ ). In particu- 
lar, the peaks in p[e) at £ = nET{n-\- 1/2), n = 0, 1, 2, . . ., 
are analogous to the De G ennes-Saint James resonances 
in J osephson junctions ( De Gennes and Saint- Jamesl 
I1963D . 

In a semiclassical description, the suppression of the 
density of states in the p-n junction can be understood 
as destructive interference of the electron-like and hole- 
like segments of a periodic orbit (solid and dashed lines in 
Fig. [77b ) . At the Fermi level, the dynamical phase shift 
accumulated in the p and n regions cancels, and what 
remains is a Berry phase shift of tt from the rotation of 
the pseudospin of a Dirac fermion. 

If the p and n regions enclose a magnetic flux as in 
the ring geometry of Fig. [28] (inset), then the Berry phase 
shift can be compensated and the suppression of the den- 
sity of states can be eliminated. The resulting flux de- 
pendence of the ground state energy E — A p{£)£ de 
(with A the joint area of the n and p regions) implies 
that a current / = dE/d^ will flow thro ugh the rin g 
in equilibrium, as in a Josephson junction ( Imrvl llQQTi ). 
According to Eq. ([5^ . the order of magnitude 



(37) 



Jo = {e/h)El/S ^ {e/ti)NEj 



of this persistent current is set by the level spacing S = 
{Apo)~^ and by the Thouless energy Et = hv/nr — NS 
in the ring geometry (of radius r and width w <^ r, sup- 
porting N — AUqw/ttTiv ^ 1 propagating modes). Be- 
cause of the macroscopic suppression of the density of 
states, this is a macroscopic current — larger by a factor 
N than the mesoscopic persistent current in a ballistic 
metal ring ( Biittiker. Imrv. and Landaueii [19831 : [lmrv[ . 
11997). 

Fig. [55] plots /(<&) for an abrupt p-n junction 
i n an A^-mode ring without intermode scattering 
( Beenakker et al.l . [20071 ). The maximal persistent cur- 
rent is Ic ~ 0.2 /q. Up to a numerical coefficient, this 
result for h is the same as the critical current of a bal- 
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